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Preface 



The aim of these lectures is to give a self-contained introduction to 
nonrelativistic potential models, to their formulation as well as to their 
possible applications. At the price of some lack of (in a mathematical 
sense) rigorous derivations, we try to give a feeling and understanding 
for the simplest conceivable method to extract the explicit form of the 
forces acting between quarks from the interplay between experimental 
observations and theoretical considerations. According to this spirit, we 
demonstrate, in detail, how to obtain the underlying Hamiltonian and 
how to determine the Lorentz structure of the quark-(anti-)quark in- 
teraction potential from well-established experimental facts. 

Vienna, December 1995 

Wolfgang Lucha 
Franz F. Schoberl 
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Chapter 1 

Nonrelativistic Potential 



In principle, the appropriate framework for the description of bound 
states within relativistic quantum field theories is the Bethe-Salpeter 
formalism. There are, however, some circumstances which are opposed 
to this. The Bethe-Salpeter equation cannot be solved in general. The 
interaction kernel entering in this equation is not derivable from QCD 
either. The propagators of the constituents have to be approximated 
by their free form, the involved masses, however, being interpreted as 
effective ("constituent") ones. So, even if one is wilhng to put up with 
the complexity of the Bethe-Salpeter formalism, it is hard to obtain 
information from this approach. 

The alternative which comes closest to one's physical intuition is the 
description of bound states with the help of the Schrodinger equation 

as a 

where the nonrelativistic Hamiltonian for a quantum system consisting 
of two particles with masses mi and m2, respectively, which interact via 
some potential V(x) is given in the center-of-momentum frame by 

= mi + 7712 + h V{yi) ; 

here /i denotes the reduced mass, 

mi m2 
/i = . 

mi + m2 

Our main task is simply to find that potential ^(x) which describes 
the interaction of the two particles constituting the bound state under 
consideration. By investigating the corresponding scattering problem of 
the involved bound-state constituents, the perturbatively accessible part 
of this potential may be derived according to the following recipe (for 
details see, for instance, Refs. H, 0]): 
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1. Compute the scattering amplitude Tfi, which is defined in terms of 
the S-matrix element introduced in Appendix 

= (f, out|i, in) , 

by the decomposition 

for the elastic scattering process i ^ f in lowest non-trivial order 
of perturbation theory, the so-called "first Born approximation." 

2. Perform the nonrelativistic limit, realized by the vanishing of the 
momenta p of the involved bound-state constituents; we indicate 
this limit rather symbolically by 

3. Obtain the configuration-space interaction potential sought after, 
y(x), as the Fourier transform of the above scattering amplitude 

V(x) = -{27Tf J d^kexp{-ik-x.)Tn{k) . 

For the sake of simplicity, we split off all the normalization factors 
of the one-particle wave functions, given for a fermion of mass m and 
kinetic energy 



by 



1 



m 



we thereby define a quantity t according to 



In the framework of this nonrelativistic treatment, the short-range 
part of the quark-antiquark potential (which is of perturbative origin!) 
will be determined from quantum chromodynamics (QCD) according 
to the above prescription. The shape of the long-range, confining part 
of the quark-antiquark potential (which is of nonperturbative origin!) 
will be obtained from the analysis of the possible Lorentz structures of 
the potential, its coordinate dependence from the comparison of some 
of the resulting predictions with experiment. 



1.1. NONRELATIVISTIC LIMIT 
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1.1 Nonrelativistic limit 

In the nonrelativistic limit, our whole formahsm, so to say, "collapses" 
to an extremely simple one: 

• The relativistically correct expression for the one-particle kinetic 
energy, 

Ep = + m2 , 

reduces to 

Ep = m . 

• The Dirac spinors u{p,a) and v{p,a) describing fermions of mass 
m, four-momentum p, and spin polarization cr. 



u{p,a) = 



v{p,a) = 



Eri + m 



m 



( 1 \ 
cr • p 



Xa 



E„-\- m 



m 



\ Ep + m J 

( \ 

Ep + m 

V 1 / 



Xa ^ Xa=-^(^2Xa 



where cr = {ai,i = 1,2,3} are the three Pauli matrices 



ai 



II 1 
1 



(72 



II -I 

i 



0-3 



I 
-1 



reduce to the nonrelativistic spinors 



Xa 





Act 



• The normalization of the Dirac spinors adopted by us, 

u\p, cr) u{p, r) = v\p, cr) v{p, t) = ^ S^t , 

m 

reduces to the nonrelativistic normalization 

u\(t)u{t) = xiXr = ^ar , 
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1.2 Static potential in quantum electrodynamics 

Let us illustrate the very simple procedure outlined above by applying 
it first to electron-positron scattering: 



e (gi,Ti) 



7 



e7i> 



e {pi,ai) ' ' e+(p2,cT2) 

Figure 1.1: Electron-positron scattering, one-photon exchange graph. 

The interaction term in the Lagrangian of quantum electrodynamics 
(QED) for the coupling of a fermion with electric charge e, described 
by the Dirac spinor field ip{x)^ to the photon field A^{x) is well known: 

jC^^^{x) = e^^J{x)J^^|J{x)A^'{x) . 

In lowest non-trivial order of the perturbative loop expansion, just 
two Feynman diagrams contribute to the scattering amplitude Tfi for 
elastic electron-positron scattering: 

• the exchange of a single photon, 7, between electron and positron, 
as depicted in Fig. [TTT|; and 

• the annihilation of the electron-positron pair into a single photon, 
7, followed by a subsequent creation of an electron-positron pair 
by this single photon, as depicted in Fig. [O . 



With the help of the Feynman rules given for an arbitrary, that is, in 
general, non-Abelian, gauge theory in Appendix |5, it's straightforward 
to find the corresponding scattering amplitude Tfi: 



1.2. STATIC POTENTIAL IN QUANTUM ELECTRODYNAMICS 




e+(g2,T-2) 



7 




Figure 1.2: Electron-positron scattering, pair annihilation graph. 

• The contribution of the one-photon exchange graph in Fig. [TTT] to 

our scattering amplitude t reads 

e2 _ 

texch = --^u{qi,Ti)'y^u{pi,ai)v{p2,(J2)Yviq2,r2) , (1.1) 

where k denotes the involved momentum transfer, 

k=pi-qi = q2-p2 ■ 
The square of this momentum transfer, 

e = {pi-qif 

which enters in the denominator of the scattering amplitude texch, 
reduces in the nonrelativistic limit to 

= -k^ . 

The spinor factors u^^ u and v'y^v may be evaluated very easily: 

— For our particular choice for the normalization of the Dirac 
spinors u and we have, in the nonrelativistic limit, 

w(ri)7ow(cri) = u\ti) u{ai) = Sr.ai , 

^((^2)70^^(7-2) = v\a2)v{T2) = 5r^a2 ■ 
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In the Dirac representation, the Dirac matrices 7^ = {7^, 7} 
are explicitly given by 



70 



7 



1 

-1 

cr 

-cr 



Inserting these explicit representations of the Dirac matrices, 
we obtain, in the nonrelativistic limit, 

u{Ti)ju{ai) = w^(ri)7o7w(cri) 



1 \ / cr 



0-1-^0 '^^^^^ 



a 



(Xn^O 




cr\ Xa, 
o- \ 



and, similarly. 



v{(72)jv{t2) = . 

Accordingly, the scattering amplitude texch of Eq. (|TTT|) reduces to 



e2 



^exch — ^'''1^1 ^T2<^2 • 

• The contribution of the pair annihilation graph in Fig. |i.2| to our 



scattering amplitude t reads 

tann = 7 ] u{qi, Ti) J^v{q2, T2) v{p2, (^2) 7^ ^{pi, (7i) . 

Here, the total momentum P of the system under consideration, 

P =Pi+P2 = qi^q2 , 

enters in the denominator of the scattering amplitude ^ann- In that 
case, however, the square of this total momentum. 
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reduces in the nonrelativistic limit to 

Thus, compared with the contribution to the scattering amplitude 
t arising from one-photon exchange, ^cxch, the contribution to the 
scattering amplitude t arising from pair annihilation, tann, will be 
of the order of 

k2 



This observation indicates that the annihilation contribution tann 
represents, in any case, already some relativistic correction to the 
exchange contribution texch, for which there will be no room at all 
within a purely nonrelativistic investigation and which, therefore, 
has to be neglected for the present discussion. 

We are unambiguously led to the conclusion that, in the nonrelativistic 
limit, only the one-photon exchange graph contributes to the T-matrix 
element for elastic electron-positron scattering: 

1 1 

According to step 3 of our procedure, the interaction potential V{x.) 
is obtained as the Fourier transform of the T-matrix element T^. Since, 
at present, we are exclusively interested in the nonrelativistic limit, we 
shall obtain in this way only the nonrelativistic (or static) part 14[r(x) 
of the potential: 

V^nr(x) = -(27r)^y' Aexp(-2k-x)rfi(A;) 
= — / (i^A;exp(— ik • x) t 
1 g2 /^3,.exp(-ik-x) 



(27r)3 J k2 

The result of the required integration may immediately be written 
down: 

1. The integral is obviously invariant under rotations. Consequently, 
it has to be some function $ of the radial coordinate r = |x| only: 



/ 



3 exp(-zk • x) 
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2. For dimensional reasons, this function $(r) has to be proportional 
to the inverse of r: ^ 

$(r) oc - . 
r 

These considerations justify the ansatz 

,3, exp(-ik-x) A 



(27r)3^ k2 r ' 

with some dimensionless constant A. We determine the constant A by 
applying the Laplacian A = V • V to both sides of this ansatz: 

• For the left-hand side, we find 

1 A /■ T-^, exp(-zk-x) 1 /• o, , , 

= -<5P)(x). 

• For the right-hand side, upon remembering the relation 

Ai = -47r(5(3)(x) , 
r 

we find 

AA- = -47rA(5(^)(x) . 
r 

By comparison, the dimensionless proportionality factor, A, is pinned 
down to the value 

A = ^. 

With due satisfaction, we realize that, by following step by step our 
general prescription given in our introductory remarks to this chapter, 
one is indeed able to recover, from the nonrelativistic limit of the Born 
approximation to the T-matrix element for (elastic) electron-positron 
scattering, the static Coulomb potential of quantum electrodynamics: 

T/QED/ X ^ 

47rr 

or, with the usually employed definition 

of the electromagnetic fine structure constant. 
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1.3 Static potential in quantum chro mo dynamics 

The overwhelming success in the case of quantum electrodynamics has 
contributed to enhance our confidence in our prescription of extracting 
the (perturbatively accessible part of an) effective interaction potential 
from the relevant elastic-scattering problem. Hence, we do not hesitate 
to apply this procedure also to the case of quantum chromodynamics. 

The relevant situation for the determination of the potential which 
describes the quark forces acting within mesons is the quark-antiquark 
scattering 

where the indices j, . . . = 1, 2, 3 denote the colour degrees of freedom 
of the involved quarks. 

According to our brief but nevertheless comprehensive — not to say, 
exhaustive — sketch of quantum chromodynamics given in Appendix 0, 
the coupling, with the interaction strength ^^s, of a quark represented 
by the Dirac spinor field qi{x), to the gluon fields a = 1, 2, . . . , 8, 

is described by the interaction Lagrangian 

A"- 

d^^^{x) = gs q^{x) 7^ qj{x) C^ix) , 

where A^, a = 1, 2, . . . , 8, are the eight Gell-Mann matrices; an explicit 
representation of these matrices may be found in Appendix |U.2| . They 
serve to construct a fundamental (three-dimensional) representation of 
the generators of the gauge group SU(3) of quantum chromodynamics: 

fund 2 

Because of the (structural) similarity of the interaction Lagrangians 
jCi of quantum electrodynamics, >Cj^^^, and quantum chromodynamics, 
>C^^^, again only two Feynman graphs potentially contribute, in lowest 
non-trivial order of the perturbative loop expansion, to the scattering 
amplitude for elastic quark-antiquark scattering: 

• the exchange of a single gluon, G^, between quark and antiquark, 
as depicted in Fig. [C3| ; and 



the "annihilation" of the quark-antiquark pair into a single gluon, 
Gq, followed by the subsequent creation of a quark-antiquark pair 
by this single gluon, as depicted in Fig. [TT^. 
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meson 




meson 



Figure 1.3: Quark-antiquark scattering, one-gluon exchange graph. 



meson 



Qk{qi,ri) 




G 




Qj{P2,cr2) 



meson 



Figure 1.4: Quark-antiquark scattering, pair annihilation graph. 
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There is absolutely no need to calculate the corresponding T-matrix 
elements Tfi once more. Comparing the above interaction Lagrangians 
jCi of quantum electrodynamics, and quantum chromodynamics, 

jC^*^^ we realize that — loosely speaking — we may obtain the transition 
amplitudes required in the present case from the ones computed in the 
previous section by simply replacing in the latter the electric charge e 
by the expression 

However, we have to take into account that, according to the famous 
confinement hypothesis, all the quarks inside a hadron form necessarily 
a colour-singlet state. Consequently, we feel obliged to amend — which, 
as there is no danger of confusion, we do without change of notations — 
the transition amplitudes Tn for the scattering of free particles by the 
(normalized) meson colour wave functions 

J_. 
V3 ' 

which means, in fact, nothing else but an appropriate average over the 
colour degrees of freedom. 

With this proviso, the T-matrix element Tfi in question is found as 
follows: 

• The contribution of the one-gluon exchange graph in Fig. \L.'6\ to 

our scattering amplitude t reads 

^exch = ^u{qi,Ti)-f^u{pi,ai)v{p2,(72)Yv{q2:r2) , 

where k denotes again the involved momentum transfer, 

k = pi - qi = q2- P2 , 

or, after the announced multiplication by those meson colour wave 
functions, 

^2 1 1 \a \a 

+ - _^_La -La 
'"'^ ~ /c2 V3 V3 2 2 

X u{qi,Ti)-i^u{pi,CFi)v{p2,(J2)l^v{q2,r2) . 

With one of the relations given in Appendix |D.3| , the colour factor 
stemming from the exchange graph yields 

13 13 8 \a \a„ 



\/3 i^j^i a/3 k,£=l a^l 2 2 



meson colour wave functions "(colour charge)^" /(/g 
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X 3 8 

— 77^ "12 Ki Kk 

16 

_ 4 
~ 3 ■ 

Hence, the one-gluon exchange contribution igxch to the scattering 
amplitude t is given by the expression 
4 _ 

which reduces in the nonrelativistic hmit to 

3 k2 



^exch — o 1 ^Tiai ^T2(J2 



• The contribution of the pair annihilation graph in Fig. |i.4| to our 
scattering amplitude t reads 

4nn = 7 -TT U{qi, Ti) -ff,v{q2, T2) v{p2, (J2) 1^ u{pi, di) , 

or, after the announced multiplication by those meson colour wave 
functions. 



^2 -1 1 \fl \a 

(1^1+1.2)2 V3 ^^'x/3''' 2 2 



X w(gi,ri)7^i;(g2,T-2)?;(i?2,cr2)7^wbbcri) . 

However, this annihilation contribution ^ann vanishes identically: 

1 3 A"- 1 



\/3 ij=i 2 2\/3 
= . 

From the physical point of view, the interpretation of this, at first 
sight slightly surprising, phenomenon is rather evident: the gluon, 
as a colour octet, that IS, as a particle which transforms according 
to the eight-dimensional adjoint representation of SU(3)c, has no 
means to couple to a colour singlet, like any bound state of quarks, 
without violating thereby the conservation of colour demanded by 
the exact invariance of quantum chromodynamics with respect to 
the colour gauge group SU(3)c. 
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Collecting all our above findings, we may state that, in lowest order 
of the perturbative loop expansion, the T-matrix element Tfi for elastic 
quark-antiquark scattering within mesons receives only a contribution 
from the one-gluon exchange graph: 

1 1 4 ^2 

Consequently, the quintessence of the present consideration is: we may 
recover the (perturbatively accessible part of the) effective interaction 
potential operative in quantum chromodynamics from its counterpart 
in the case of quantum electrodynamics by simply replacing the square 
of the electric charge e by the factor 

4 2 

Hence, the short-distance part of the static quark-antiquark potential, 
arising from one-gluon exchange within mesons, is of Coulombic shape: 

3 47rr 

or, with the usually employed definition 

a =i. 



of the strong fine structure constant. 



6 r 



1.4 Lorentz structure of an interquark interaction 

At this stage, in order to seize hold of the nonperturbative contribution 
to any effective potential, we embark on a rather general investigation. 

Quite generally, the T-matrix element Tfi for the elastic scattering of 
some generic fermion and the corresponding antifermion both of 
them of mass m, 

is, apart from the overall normalization factors of the one-particle wave 
functions, which we always split off by the definition 

_ 1 

^' (2^ V^pi ^.1 ^ 



of our scattering amplitude the product of 
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• two bilinears of Dirac spinors of the form u{qi, ri) Fi u{pi, ai) and 
v{p2, (^2) ^2v{(i2^ T2), where Fi and F2 represent some (unspecified) 
Dirac matrices, and 

• some interaction kernel K which, a priori, may depend on all four 
external momenta Pi, P2, Qi, Q2, 

K = K{pi,p2,qi,q2) , 

only subject to the momentum conservation 

Pi + P2 = + , 

as expressed by that overall 6 function multiplying this T-matrix 
element Tfi in the standard decomposition of the S-matrix element 

Consequently, the most general ansatz for our scattering amplitude t 
reads 

t = u{qi, n) Fi u{pi, CTi) v{p2, (72) F2 v{q2, T2) K , 

with K depending on any three independent linear combinations built 
from the external momenta out of the set {pi,P2, Qi, ^2}- 

We shall constrain the T-matrix element Tfi under consideration by 
the following two, very reasonable assumptions: 

1. The T-matrix element Tfi is invariant with respect to the full set 
of transformations forming the (homogeneous) Lorentz group, that 
is, invariant under 

• proper orthochronous Lorentz transformations, 

• space reflection ("parity operation"), and 

• time reversal. 

2. The interaction kernel K entering in the T-matrix element Tfi is a 
function of only the square of the involved momentum transfer 

k=pi-qi = q2-p2 ; 

that is, 

K = K{k^) , 
which reduces in the nonrelativistic limit to 

K = K{-k^) . 
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The most general form of the scattering amphtude t consistent with 
the requirements of the above assumptions is (see, e.g., Refs. [g, g, [7j) 

t= E t^, (1-2) 

S=S,P,V,A,T 

where any particular contribution t^, is of the form 

tj: = u{qi, Ti) u{pi, ai) v{p2, (12) Te v{q2, T2) K^{k^) 
and the sum extends over the five possible Lorentz structures 
• scalar (S), 



= 1 1 , 

Ts (g) = 75 (g) 7^ , 
(g) = 7,, g) 7^ , 

rsg)rs = 7^75 7^7^ , 



• pseudoscalar (P). 

• vector (V), 

• axial vector (A), 

and 

• tensor (T), 

Here, we had to introduce the Dirac matrices 

^/ — . ^ 

= -?7°7^7^7^ 
= ^70 717273 , 

with the totally antisymmetric Levi-Civita symbol in four dimensions 

unambiguously fixed by demanding 

eoi23 = 1 , 

as well as 
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In the Dirac representation, these Dirac matrices are explicitly given by 

75 = - 



1 

1 



and 



a^' = i\ " 7 1 , z = 1,2,3 



a 



a, 



5 "5 " ? 



c^'^ = e^jk ^^j , ij,k = 1,2,3 . 

In the two preceding sections, that fermion-antifermion interaction 
was basically mediated by the exchange of some vector boson, namely, 

• the photon in the case of quantum electrodynamics or 

• the gluon in the case of quantum chromodynamics. 

As a consequence of this, in both of these cases the effective interaction 
was solely of vector Lorentz structure, 



Ti (g) r2 = ® J 



or, in other words, in the decomposition ( |1.2|) the scattering amplitude 



t received exclusively a vector contribution. The interaction kernel K 
was given by 



reducing to 



where the effective coupling strength k, stands 
• in the case of quantum electrodynamics for 



2 

K = e , 

• in the case of quantum chromodynamics for 

4 2 

Fourier transformation then resulted in the static Coulomb potential 



47r r 
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Needless to say, in general the various contributions entering in the 
decomposition ([TT2j) of our scattering amplitude t will not arise from 



the exchange of a single particle representing some fundamental degree 
of freedom of the underlying quantum field theory. Rather, these terms 
have to be interpreted as due to only an effective exchange of a particle 
of the appropriate behaviour under Lorentz transformations. 

We should be prepared to the fact that — in the course of evaluating 
below the various terms t^, contributing to the T-matrix element Tfi — 
we shall encounter expectation values of the Pauli matrices of the form 
Xn^ Xcri Q-iid Xct2^ ^ Xt2 - shall cast the second of these expressions, 
which involves two-component spinors Xa representing the spin degrees 
of freedom of antifermions, defined by 

Xl = X*a , 

with the help of the identity 

T 

(72<T(72 = —Cr 

into an equivalent form which involves only the two-component spinors 
Xa pertaining to fermions: 



Xa <^Xt = [Xaf^X 



T 



Xa (72(^CF2X 



-Xt crXa 



We shall find it convenient to abbreviate the expectation values of the 
Pauli matrices cr by introducing the shorthand notation 

CT, = -X^>X^, (1.3) 

Xt2 ^ Xoi ■ 

Now, what we really have to do when trying to follow the steps given 
in our prescription for the derivation of that (perturbatively accessible 
part of an) effective interaction potential from the underlying quantum 
field theory by considering the relevant elastic-scattering problem may 
be phrased in the following way: 
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1. Calculate the expectation values of the considered S operator, S'fi, 
or T operator, Tfi, with respect to those Fock-space states usually 
employed in quantum field theory. 

2. Look upon these matrix elements as the expectation values of that 
interaction potential you are searching for, V(x), with respect to 
the quantum-theoretical bound states and extract this interaction 
potential by "factorizing off" all remnants of these bound states. 

In this and only this (!) sense one may adhere, when switching from the 
scattering amplitude to the interaction potential, to the identifications 
of the spin operators Si and S2 of fermion and antifermion with half of 
the expectation values <ti and 0-2 of the Pauli matrices a, respectively: 

^1 = ^ ? 



S2 "=' 

and therefore 



2 
2 



<7i • <72 "=" 4 Si • S2 . 

Adopting an admittedly rather symbolical notation, we suppress in 
the following any reference to both the momenta and the spin degrees 
of freedom of the involved particles. Nevertheless, at every moment it 
should be clear what's going on. With the above at one's disposal, the 
nonrelativistic potentials ^^^^(x), S = S, P, V, A, T, are easily found: 

Scalar: For the scalar Lorentz structure, i.e., 

Fe Fs = 1 1 , 

we find for our scattering amplitude t in the nonrelativistic limit 

= -Ksi-k') . 

Upon Fourier transformation, the corresponding static interaction 
potential ^nr(x) reads 

Vi2i^) = Vs(r) , 
where Vs{r) is defined by 

Vsir) = J d'keM-i^- x) i^s(-k') • 

Accordingly, the scalar Lorentz structure yields a "pure potential" 
term. 
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Pseudoscalar: For the pseudoscalar Lorentz structure, i.e., 

Te (g) Te = 75 7^ , 
we find for our scattering amplitude t in the nonrelativistic limit 

= 7o 75 w "U^ 7o 7^ "y i^p(A;^) 
= . 



Consequently, in the nonrelativistic limit, the contribution of this 
pseudoscalar Lorentz structure vanishes: 

<^(x) = 0. 
Vector: For the vector Lorentz structure, i.e., 

Fs Fs = 7a. t'' , 

the nonrelativistic limit of our scattering amplitude t has, in fact, 
already been calculated in Section [TT2|: 

Upon Fourier transformation, the corresponding static interaction 
potential V^^'*(x) reads 

where Vy{r) is defined by 



Vyir) = 



1 r 

— / d^k exp(-i k • x) Kyi-k' 



Accordingly, the vector Lorentz structure yields, very similarly to 
its scalar counterpart, a "pure potential" term. 

Axial vector: For the axial vector Lorentz structure, i.e., 

Fe(^Fe = 7/i 75(^)7^7^ , 

we find 
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• for // = 



wyo u vj^ V = 75 w 7^ v 



X 




= 

and 
• for // = i 

U'J ■ V 'J V = ix^ 70 7 75 u • -y^ 7^ 7 7^ -u 



= (x^o)( 



W X 

-(T ) [0 

-a W 

-^^ j 



= — CTi • (72 



and therefore for our scattering amplitude t in the nonrelativistic 
hmit 

= (Ti • <72i^A(-k^) • 

Upon Fourier transformation, the corresponding static interaction 
potential V^^''(x) reads 

y,J^)(x)=4Si-S2yA(r) , 
where Vp^{r) is defined by 

VA{r) = J d'k exp(-z k • x) KA{-k') . 

Accordingly, the axial vector Lorentz structure entails an effective 
spin-spin interaction. 
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Tensor: For the tensor Lorentz structure, i.e., 

Te (K) = ^ a^^ (8) (7^" , 

we find 

• for ij, = 0,1/ = i 



= 



aA 



(J, 



and 

• for fi = i, u = j 



= e^jfc(x^O) 



CTfc 



X 



-CTk \0 



X 



^iji (0, 







-^tW I t 



= -2y^ 



X CTkXX"^ CFk X^ 

= -'^^ (^x-x^'^ ^X^ 

and therefore for our scattering ampUtude i in the nonrelativistic 
hmit 

1 



1 



u aoi u V a""' V -\ — u aij uva^^ v 
2 



= CTi • (72 i^T(-k^) 



Upon Fourier transformation, the corresponding static interaction 
potential V^r (x) reads 
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where Vr^r) is defined by 

Vrir) = J d'k exp{-i k • x) KT{-k^) . 

Accordingly, the tensor Lorentz structure entails also an effective 
spin-spin interaction. 

Table summarizes our findings for the contributions of the various 
possible Lorentz structures to the effective interaction potential in the 
nonrelativistic limit. 

Table 1.1: Nonrelativistic interaction potential V^r arising effectively from the various 
conceivable Lorentz structures Fs^Ts of an arbitrary fermion-antifermion interaction 



Lorentz structure Fs <S) Fs static potential 

scalar 1 (8) 1 ^s(^) 

pseudoscalar 75 ® 7^ 

vector 7^ (g) 7^^ W(^) 

axial vector 7^ 75 (g) 7^^ 7^ 4 Si • S2 VA(r) 

tensor | cr^i. cr^'' 4Si-S2VT(r) 



The total spin S of the respective bound state under consideration 
is clearly given by the sum of the spins Si and S2 of its constituents: 

S = Si + S2 . 

Upon squaring this relation, 

S^ = (Si + S2)2 

= S^ + S2 + 2 Si • S2 , 

we may express the product Si • S2 of the two spins Si and S2 in terms 
of the squares of Si, S2, and S, 

Si-S. = i(S^-S?-S^) , 

and, therefore, its expectation values by the corresponding expectation 
values of S^, S2, and S^: 
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Accordingly, expressed in terms of the quantum numbers 5, 5*1, and 5*2 
of the spins S, Si, and S2, respectively, the expectation values (Si • S2) 
of the product of the spins Si, S2 of the bound-state constituents read 

(Si -82) = ^ [S (5 + 1) - Si {Si + 1)-S2 {S2 + 1)] . 
For fermionic constituents with spin 

^1 = ^2 = ^ , 



we have 



and therefore 



Si (Si + 1) = S2 (52 + 1) = I 

(Si.S2) = i5(5 + l)-?. 
Moreover, for fermionic constituents with spin 

^i = 5*2 = - , 

the quantum number S of the total spin S may accept precisely either 
of two values: 

• S = 0, which corresponds to some spin singlet, hke the pion or the 
77 meson in the case of light quarks, or the r]c in the charmonium 
system. 

• = 1, which corresponds to some spin triplet, like the cj, and (f) 
mesons in the case of light quarks, or the J/ip in the charmonium 
system, or the T in the bottomonium system. 

This implies for the eigenvalues S {S + 1) of the square S^ of the total 
spin S: 



5(5+1) = < 



for spin singlets, i.e., 5 = 0, 
2 for spin triplets, i.e., 5 = 1. 

Accordingly, the expectation values (Si • S2) of the product of the spins 
Si, S2 of the bound-state constituents are finally given by 

3 



(Si • S, 



for spin singlets, i.e., 5 = 0, 

* (1-4) 

+- for spin triplets, i.e., 5=1. 
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The first and simultaneously most important lesson to be learned from 
the above is that, for these two possible values of the quantum number 
S of the total spin S of the bound state, the spin-spin interaction term 
will contribute necessarily with opposite signs. 

Collecting all previous results, the following picture emerges for the 
nonrelativistic limit of the most general effective fermion-antifermion 
interaction potential T4^r(x): 

• Both scalar and vector Lorentz structures lead to "pure potential" 
terms: 

^NR (x) = for S = S,V. 

• The contribution of the pseudoscalar Lorentz structure vanishes:^ 

^iS(x)=0. 

This circumstance provides, for instance, a very compelling reason 
for the (relatively) weak binding of deuterium in nuclear physics: 
The interaction between nucleons, that is, protons and neutrons, 
is generally accepted to be dominated by one-pion exchange. Since 
the TT meson is a pseudoscalar meson, the Lorentz structure of its 
coupling to the nucleons has to be also of pseudoscalar nature in 
order to form an interaction Lagrangian which is a Lorentz scalar. 
This fact implies that only the relativistic corrections arising from 
one-pion exchange can be responsible for the binding of a proton 
and a neutron to the deuterium. 

• Both axial vector and tensor Lorentz structures contribute merely 
to the spin-spin interaction term: 

^NR(x) = 4Si-S2VE(r) for S = A, T . 

As a consequence of Eq. ([TT^), in this case we will obtain a binding 
force between fermion and antifermion only for either of the above 
two possible values, S = and 5 = 1, of the quantum number S 
of the total spin S of the bound state. 

Hoping that the empirically observed hadron spectrum will provide 

some restrictions on the allowed effective quark-antiquark interaction, 

we now confront the above picture of general findings with experiment: 

^ In the next — i.e., first non-trivial — order of tlie present nonrelativistic expansion, 
tliis pseudoscalar Lorentz structure contributes to tlie spin-spin and tensor interaction 
terms. 
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• Already the mere existence of strongly bound mesons forbids the 
pseudoscalar Lorentz structure to play any significant role within 
some phenomenologically acceptable quark-antiquark interaction 
potential. 

• The existence of both pseudoscalar mesons (like tt, 77, and 77') and 
vector mesons (like p, u, and (f)), all of which are bound states of 
a quark-antiquark pair with vanishing orbital angular momentum, 
implies that the actual quark-antiquark forces must be described 
by an interaction potential which yields binding for 5* = as well 
as 5 = 1. Obviously, this fact rules out both the axial vector and 
tensor Lorentz structures as the predominant contribution to any 
realistic quark-antiquark interaction potential. 

In other words, the theoretically predicted particle spectra would look 
very different from the experimentally measured ones if the dominant 
terms in the effective quark-antiquark interaction potential would not 
be just some linear combination of vector and scalar Lorentz structure. 
Therefore, our conclusion has to be: 

The Lorentz structure of the quark-antiquark interaction 
is dominated by the vector 7^ (8) 7^ and / or the scalar 1 (8) 1 , 
both of which lead in the nonrelativistic limit to so-called 
pure potential terms. Thus the static interaction potential 
^nr(^) must be the sum of merely the contributions of the 
vector — Vv{r) — and the scalar — Vs(r) — Lorentz structures: 

V^R{r) = Vyir)+Vsir) . 



Chapter 2 

Relativistic Corrections 



Beyond doubt, the next logical step must be to improve the up-to-now 
entirely nonrelativistic formalism by taking into account all relativistic 
corrections. In principle, one encounters no particular difficulties when 
trying to take into account (at least, at some formal level) the complete 
set of relativistic corrections to the effective interaction potential |5j] . 

For the moment, however, we intend to be somewhat more modest, 
and this even in two respects: 

1. We shall calculate these relativistic corrections only up to second 
order in the absolute value = |v| of the generic relative velocities 

V = — 
Ep 

of the bound-state constituents, that is, only up to order 

2 

P 

which, since up to this order the relativistic kinetic energy 



Ep = yp"^ + m'^ 
may be approximated at this place by 

is equivalent to 

2 y 



We shall consider only the spin- dependent contributions to these 
relativistic corrections. These spin-dependent interactions control 
the fine and hyperfine level splittings of the bound-state spectra. 
The spin-independent interactions may be obtained, with slightly 
more effort, along similar lines [0. 
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2.1 Spin-dependent corrections 



We shall be interested in all the spin-dependent relativistic corrections 
to the static interaction potential Vnr(x) up to order 

^2 



Therefore, we must focus our attention to those terms in the transition 
amplitude Tfi which involve expectation values of the Pauli matrices a. 
One may simplify this task considerably by the following observation: 
we are entitled to approximate the relativistic kinetic energy 



by the lowest-order term 



Ep^m 



of its nonrelativistic expansion in two places, namely, 
1. in the Dirac spinors u{p^ a) and v{p., cr), 

/ 1 \ 



u{p,a) = 



£■-0 + 777. 



m 



Ej) + m 



m 



• P 

\ Ep + m J 

( \ 

Ep + m 

1 



V 



C 



-^(^2X1 



I 



which then (only for the purpose of the present analysis!) become 

/ 1 \ 



• p 

\ 2m J 
( ^ ^ 

2 777 

V 1 J 



Xa 



(2.1) 



(7 ' 



2. in that general relationship between the T-matrix element Tfi and 
our scattering amplitude t, 



Tfi=: 



1 



m 



(27r)6 ^Ep^ Ep^ E,^ E,, 



t 



which then becomes 



Tfi = 



(27r)^ 
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In the course of calculating the scattering amplitude we may take 
advantage of two trivial simplifications: 

1. In order to get the interaction potential ^(x), we have to consider 
the scattering amplitude t only in the center-of-momentum frame, 
defined by the vanishing of the total momenta Pj and Pf of initial 
state i and final state f, respectively: 

Pi = Pi + P2 = qi + q2 = Pf = . 

Consequently, this scattering amplitude t will depend only on the 
involved momentum transfer 

k = pi - qi = q2 - P2 
and on the relative momentum 

P = Pi = -P2 

2. The Pauh matrices a = {dj, i = 1,2, 3} satisfy both 

• the commutation relations 

WtjCTj] = 2ieijk(7k 

and 

• the anticommutation relations 

{(Ti,(7j} = 2 5ij . 

Adding up these two relations, the product Oi aj of any two Pauli 
matrices is given by 

= 6ij + icijkdk ■ 

By application of this relation, any product of (two or more) Pauli 
matrices may be reduced to an expression which involves no more 
than at most one Pauli matrix. 

Performing the Fourier transformation as demanded by that general 
prescription briefly sketched in our introductory remarks to Chapter 
the effective interaction potential V{'x) is derived from t according to 

y(x) = -{27rf J d^kexp{-ik-^)Tnik) 
= — - — (fkexpi—ik ■ -xjt . 
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As the central result of the intended inclusion of all spin-dependent 
relativistic corrections up to order 



we shall finally end up with the generalization to arbitrary interaction 
potentials of the well-known Breit-Fermi Hamiltonian, of the standard 
form 

H = 1711 + 1712 + h V(x) , 

with fi the reduced mass of the particular two-particle quantum system 
under consideration, 

mi 777.2 

/i = . 

7771 + 7772 

Here, the interaction potential ^(x) will encompass, in addition to the 
nonrelativistic contribution Vnr(x), also all spin-dependent relativistic 
corrections Vspm(x): 

y(x) = yNR(x) + l^spin(x) . 

The set of spin-dependent relativistic corrections Vspin(x) will turn out 
to consist, in general, of some 

• spin-orbit interaction term Hj^s^ 

• spin-spin interaction term i^ss, and 

• tensor interaction term i^x; 

that is, this spin-dependent part Vspin of the interaction potential V(x) 
will read 

^spin = ^LS + ^ss + . 

Bearing in mind the outcome of our analysis of the possible Lorentz 
structure of the effective interaction in a quark-antiquark bound state 
as performed in Section |i.4| , we will treat below only the case of vector 
and scalar Lorentz structure of the static interaction potential Vnr(7"): 

VMr) = Vy{r) + Vs{r) . 

Furthermore, for the sake of simplicity, we will present in the following 
all the necessary derivations in detail only for the special case of equal 
masses of the bound-state constituents: 



777i = 7772 = ^ 
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2.2 Interaction with vector Lorentz structure 



In the case of an interaction with vector Lorentz structure, that is, for 

Ts = 7^ 7^ , 

the scattering amplitude assumes the form 

tv = ri) 7^ u{pi, CTi) v{p2, (72) 7^ v{q2, T2) KY{k^) • 

Upon inserting the Dirac spinors in Dirac representation, Eq. 
and suppressing, for the sake of notational simplicity, any reference to 
the spin degrees of freedom, we find for this scattering amplitude 

tv = u{qi)-i^u{pi)v{p2)l^v{q2)KY{k^) 
= [u\qi) u{pi) v\p2) v{q2) 



■ qi 

2 m , 



( ^ \ 

(T ■ Pi 

V 2m / 



cr • q2 



2 m 

V 1 ) 



X 



X 



1 1. 



fT ■ qi 

2 m , 



X 



o- 

a 

cr 
o- 



/ 1 \ 



cr ■ pi 



V 2 



X 



m 



J 



/ o- • q2 \ 

2 m 

V 1 / 



X 



V 2m ' , 
X + T— 9 X^ (cr • qi) (cr • pi) x 
2 X'^ (cr • P2) (cr ■ q2) X 



X 



X^^ X' + 



4 m^ 
1 



4m^ 



1 



4m^ 



o- (o- • pi) X + X^ (^^ • qi) cr X 



• x'^ {cT ■ P2) o- + X^^ (^ ■ q2) x1 1 i^v(-k2) 

= |l + {X^ (^^ • qi) (cr • Pi) X + X'^ {o- ■ P2) (cr • q2) 

t 4m^ ^ 

- x^[cr(cr-pi) + (o--qi)cr]x 

• x'^ [(tr • P2) + (tr • q2)] xl} i^v(-k2) • 
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Dropping all contributions to spin-independent relativistic corrections 
at the very instant they show up, and recalling the abbreviations 
this scattering amplitude simplifies, in the center-of-momentum frame, 
to 

+ [{2p - k)i ie^jkkjGik] [{2p - k)^ + i Ctim h a2m]}'^ i^v(-k^) 



1 + [3 1 e^jk Pt kj {aik + a2k] 



CTi ■ CT2) + (CTI • k) (0-2 • k)] ) K 



Fourier transformation appears to be the appropriate tool to obtain 
the resulting effective interaction potential: 

• Already from the very beginning, all interaction kernels Ky, have 
been (implicitly) assumed to depend, in the nonrelativistic limit, 
only on the modulus |k| of the momentum transfer k. Therefore, 
all the corresponding static interaction potentials V^r (x) have to 
be spherically symmetric: 

VY{r) = — / Sk exp(— i k • x) Ky 



(27r) 

• Denoting the first and second derivatives of any static interaction 
potential with respect to the radial coordinate r = |x| by primes, 
one finds 



(27r 



^ld^kkjexp{-ik--K)Ky{-k^) = -iVjVy{r) 



r 



/ d^k k^ exp(-z k • x) Ky{-\^) = AVyir) , 
and, with the help of an identity proven in Appendix 0, 
— g / d^k ki kj exp{—i k • x) i^v(— k^) 



V,V,W(r) 

1 



Xi Xj 1 

~ 3 
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Consequently, the spin-dependent relativistic corrections V^^^°^ for the 



case of a vector Lorentz structure of the effective fermion-antifermion 
interaction become 
3 



vector 



* spin 



+ 
+ 



2 r 
2 



(xxp)-S^(r) 



3 m? 
1 



Si-S2AVv(r) 
(Si • x) (S2 • x) 1 ■ 

o 77 • 02 



lvUr)-V^{r) 
r 



Herein, it is straightforward to identify, in full accordance with the 
previously announced decomposition of the spin-dependent relativistic 
corrections ygpin, when specified to the case of vector Lorentz structure, 

T /vector rrvcctor , rrvector , rjvector 
V/_;„ = i^Ls + ^SS + 5 



spin 

• the spin-orbit term 



rrvcctor 



2m^ r 
3 



(xxp) -814(0 
L -814(0 , 



2m^ r 

with the relative orbital angular momentum 

L = X X p 
of the bound-state constituents; 
• the spin-spin term 

2 



rrvector 
^SS 



3 m" 



8i-S2Ayv(r) ; 



and 

• the tensor term 



rrvector 



1 



Si • x) (82 • X 



Si -8, 



1 



s 



12 



r 



I2m? 

with the shorthand notation 



Su = 12 



(Si • x) (82 • x) 1 ■ 

o 77 '-'I ■ "^2 



for the spin-dependent factor. 
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2.3 Interaction with scalar Lorentz structure 



In the case of an interaction with scalar Lorentz structure, that is, for 

Te ® Ts = 1 (8) 1 , 
the scattering amphtude { \r?2\ ) assumes the form 

Upon inserting the Dirac spinors in Dirac representation, Eq. 
and suppressing, for the sake of notational simplicity, any reference to 
the spin degrees of freedom, we find for this scattering amplitude 

= u{Qi)u{pi)v{p2)v{q2) Ks{k'^) 
= u^{qi) 70 u{pi) v\p2) 70 ^^(^2) i^s(fc^) 

/ 1 \ 



2m , 



1 






-1 



1 
Pi 

V 2 



X 



m 



J 



X X 



V 2m ' , 



1 






-1 



/ (T ■ q2 \ 



2 m 
1 



X'K, 



X^ X 



4 m" 



X 



X^ (cr • qi) (cr ■ pi) x 



2 X'^ {f^ ■ P2) (cr • q2) 



4 m 



4 



X^ (cr • qi) {cr • pi) x + x'^ (cr • P2) (cr • ^2) X^ 



X A's(-k' 

Dropping all contributions to spin-independent relativistic corrections 
at the very instant they show up, and recalling the abbreviations 
this scattering amplitude simplifies, in the center-of-momentum frame, 
to 

= - - -r-^{^(^tjkPikjaik + ie,jkPikja2k) \ 



1 



4 m^ 
1 

4 m^ 



i ^ijk Pi kj (cTifc + a2k) 



Once more, Fourier transformation appears to be the adequate tool 
to obtain the resulting effective interaction potential: 
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• For the same reason as before, the static interaction potential has 
to be again spherically symmetric: 



1 r 



• The corresponding spin-dependent relativistic corrections may be 
found with the help of 

^ld^kkjexp{-ik--K)Ks{-k^) = iVjVsir) 



(27r)^ 



r 

Consequently, the spin-dependent relativistic corrections V^spin^"^ ^^r the 
case of a scalar Lorentz structure of the effective fermion-antifermion 
interaction become 



T/scalar rjscalar 



with the spin-orbit term 



i7£fi^^ = -— i-L-SK'(r) , 

where L denotes, as before, the relative orbital angular momentum of 
the bound-state constituents, 

L = X X p . 

Accordingly, an interaction with scalar spin structure contributes only 
to the spin-orbit term i^Ls- However, apart from possible differences of 
the two nonrelativistic potentials Vv(r) and Vs(^) of vector and scalar 
spin structure, respectively, the spin-dependent relativistic corrections 
for scalar spin structure, Hf;^'^^^, contribute with a sign opposite to that 
of the corresponding spin-dependent relativistic corrections for vector 
spin structure, Hlf^°^. Hence, assuming identical static potentials, i.e., 

Vv{r) = Vs{r) , 

the spin-orbit term Hi^^^°^, resulting from a vector spin structure, may 
be partially compensated by the spin-orbit term i7£g^^^^, resulting from 
a scalar spin structure: 

Hlt^'=-lHir°' for Vy{r)=Vs{r). 

Before trying, in Chapter 0, to write down a (physically meaningful) 
quark-antiquark interaction potential, it is advisable to "condense" all 
these results on the effective fermion-antifermion interaction potential 
to what we would like to call a "generalized Breit-Fermi Hamiltonian." 
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2.4 Generalized Breit— Fermi Hamiltonian 



The Hamiltonian containing all spin-dependent relativistic corrections 
up to order 



V 



2 2 
2_ P' , P 



is called the generalized Breit-Fermi Hamiltonian: 




[2.2) 



where 



nil '^2 
nil + ^2 



m 



= — for nii = ni2 = m 

is the reduced mass and — according to the analysis of Section [T]J — the 
static potential V^^{r) consists of a vector and a scalar contribution, 



VNR{r) = VV{r) + Vs(r) 



(2.3) 



The corresponding spin-dependent relativistic corrections read (for the 
case of equal masses mi = ni2 = m): 



• spin-orbit term: 




[2A) 



where 

S = Si + S2 

is the total spin of the bound state and 

L = X X p 

is the relative orbital angular momentum of its constituents; 
• spin-spin term: 



H. 



ss 



3 ni' 



Si-SsAFvfr: 



;2.5) 
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• tensor term: 



H' 



1 



T 



Si • x) (S2 • x) 1 



Si-S, 



Id (f 

r dr ^ dr"^ ^ 



or, with the abbreviation 
5*12 = 12 



Si • x) (S2 • x) 1 ■ 

n ^1 ■ ^2 



:2.6) 



of the spin-dependent factor, which is sometimes called the tensor 
operator. 




:2.7) 



The corresponding expressions for the general case of unequal masses 
of the bound-state constituents mi ^ 1712 may be obtained in a similar 
manner; they are collected in Appendix |F]. 

The total angular momentum J of the respective bound state under 
consideration — which constitutes, of course, nothing else but the spin 
of the corresponding composite particle — is clearly given by the sum of 

• the relative orbital angular momentum 



L = X X p 



and 

• the total spin 
of its constituents: 



S = Si + S2 
J=L+S . 



Upon squaring this relation, 

J2 = (L + S)2 

= + S^ + 2 L • S , 

we may express the product L • S of relative orbital angular momentum 
L and total spin S in terms of the squares of L, S, and J, 

L • S = i (J2 - l2 - S2) , 
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and, therefore, its expectation values by the corresponding expectation 
values of L^, S^, and J^: 

(L . S) = ^ ((J^> - (L^> - (S^» . 

Accordingly, expressed in terms of the quantum numbers 5, and j 
of the relative orbital angular momentum L, the total spin S, and the 
total angular momentum J, respectively, denoting the bound state, the 
expectation values (L • S) of the product of orbital angular momentum 
L and total spin S, originating from the spin-orbit term i^LS, Eq. 
read 



Evidently, the expectation values of the spin-orbit term i^Ls, Eq. ( ^.4|) , 
vanish for either f = 0, 

(L • S) = for i = , 

01 S = 0, 

(L • S) = for 5 = 0, 

contributing thus only for £ 7^ and ^ = 1. The above relation yields 
the explicit nonvanishing matrix elements (L • S) listed in Table \Z.i\ . 



Table 2.1: Nonvanishing spin-orbit couplings for i ^ and S 



J 


(L-S) 


i + 1 


i 


i 


-1 


£-1 


-(^+1) 



For fermionic bound-state constituents of spin 

Si = S2 = - , 

the expectation values (Si • S2) of the product of their spins Si and S2 
in the spin-spin term i^ss, Eq. ( [2.b|) , have been determined already in 
Section [TTJ: 

3 

— - for spin singlets, i.e., 5 = 0, 



(Si • S^ 



+- for spin triplets, i.e., 5=1 
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Likewise, for fermionic bound-state constituents of spin 

Si = S2 = - , 



the spin-dependent factor 5*12 in the tensor term i^T, Eq. ( [k!.V|) , may be 
rewritten as 



Si2 = 2 



2.8) 



Accordingly, the expectation values of the tensor term i^T, Eq. ( P77|) . 
also vanish for either 5 = 0, 



(5: 



12/ 



for 5 = 0, 



or £ = 0, 

(512) = for ^ = , 

the (more or less obvious) reason for the latter being the fact that the 
rotational symmetry realized in the case of £ = leads to 



which, in turn, results in a mutual cancellation of the two terms on the 
right-hand side of Eq. ( [k!.8|) . After a lengthy calculation, the following 
expression for the diagonal matrix elements of <5'i2 may be found: 



(^12) = 



;2£ + 3) (2^-1) 



(L.S)-3((L.S))^ 



which, again only for ^ ^ and ^ = 1, yields the explicit nonvanishing 
expectation values listed in Table PT^ . 



Table 2.2: Nonvanishing diagonal matrix elements of 5*12 for £ 7^ and S = 1 



J 



^+1 



12/ 



2i 



2i + 3 
2 

2(^+1) 
2i-l 



Chapter 3 
The Prototype 



We are now in a position where we may start to think seriously about 
the question of how a reahstic, that is, phenomenologicaUy acceptable, 
potential describing the forces acting between quarks might look like. 



3.1 Funnel potential 

To begin with, let's summarize our knowledge gained so far. According 
to the analysis of Section [TT^, the quark-antiquark potential ^R(r) is 



of vector and/or scalar type, 

^NR(r) = W(r) + Vs(r) . 

For short distances, the potential — arising from one-gluon exchange — is 
(essentially) Coulomb-hke, 

6 r 

For large distances, there has to exist a contribution Fconf(^) in order 
to describe colour confinement, 

K;onf(^) = ar"^ with n > , 



implying that for large r the binding force K must not decrease faster 
than 1/r: 

K = — —^coniW) = —'nar'^~^ = —n 



dr 

From the mesonic mass spectrum, the exponent n is in the vicinity of 

n ~ 1 . 

For instance, n = 2 corresponds to the harmonic oscillator and would 
lead to equidistant level spacings. Moreover, lattice gauge theories also 
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CHAPTER 3. THE PROTOTYPE 



find that VNR(r) is roughly proportional to r for large r. Consequently, 
a linear potential, 

K;onf(r) = ar , 

is beyond doubt a sensible choice for V^onf(^)- 

The funnel (or "Coulomb-plus-linear") potential (Fig. ^TT 



4 

3 r 



one-gluon exchange 



confinement 



(s.r 



fixed in this way has been the first proposed model [§| , which in spite of 
its simplicity is able to reproduce quite well the charmonium spectrum. 
In a strict sense, the momentum (-transfer) dependence 

as = as(Q^) 

of the strong fine structure constant as has to be taken into account, 
modifying thereby the Coulomb-like behaviour of the first term on the 
right-hand side of Eq. (jXTJ) . 



V{r) 




Figure 3.1: Funnel potential. 



3.2. LORENTZ STRUCTURE OF THE FUNNEL POTENTIAL 
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3.2 Lorentz structure of the funnel potential 

In order to decide whether the Lorentz structure of the funnel potential 
is a pure vector, or a pure scalar, or a mixing of both, we consider 
the P-wave spin splittings of charmonium and bottomonium, that is, 
the ratio of mass differences 



P = 



M(3P2) - M(3P] 
M(3Pi) - M(3Pc 



(3.2) 



(Recall the usual spectroscopic notation, which designates a state with 
orbital angular momentum spin 5*, and total angular momentum j 
by "^^^^Lj, where the capital L = S,P,D,F, . . . represents the orbital 
angular momentum £ = 0, 1, 2, 3, ... , respectively.) 



Table 3.1: Masses (in GeV) and the ratio p [Eq. (pj)] for the (£ = 1,5 = 1) states of 
charmonium and bottomonium [|T^ 



Level 


(cc) 


(bb) 


(bb)' 




3.4151 


9.8598 


10.2321 


^Pi 


3.5105 


9.8919 


10.2552 


^P2 


3.5562 


9.9132 


10.2685 


P 


0.478 


0.664 


0.576 



From Table |3.1| , the experimental average for p 

p = 0.48 for (cc) , 
p = 0.66 for (bb) , 
p = 0.58 for (bb)' , 

is 

Pexp ^ 0.6 . 



With the help of the generalized Breit-Fermi Hamiltonian ( |2.2| ), we 



calculate this ratio p for the potential ( |3.1| ) perturbatively. Since the 
spin-spin interaction, iifss? does not depend on the total angular mo- 
mentum, its contribution cancels in a perturbative evaluation of p. Ac- 
cordingly, p is determined by the contributions of spin-orbit term, i^Ls, 
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CHAPTER 3. THE PROTOTYPE 



and tensor term, i^T, only 
P 





^LS + 


^P2) - fPl 


^LS + 


^Pl) 






^Pi) - (3Po 


//ls + ^tPPo) 



From Tables ^7T] and |2T^, we find for the expectation values (L • S) and 
(^12): 



and 



{S 



12 



-2 


for 


^Po 


-1 


for 


^Pi 


1 


for 


^P2 


-4 


for 


^Po 


2 


for 


^Pl 


2 


for 


^P2 



• For a pure vector, i.e. 



W = ^NR , Vs = 



one obtains 



I 8 a. (r-^) +1 air 



-IN 



P 



5 2 ac (r ^) + a (r 



,-1 



which implies the bounds 



4 7 

- < p < - 



corresponding to a = and = 0, respectively, in clear conflict 
with the experimental finding pexp — 

0.6. 



• A pure scalar, i.e. 



leads to 

which also is not tolerable from an experimental point of view. 
• A vector/scalar mixing, i.e., 

Vsj = Vexch = -| — , Vs = Vconf = a T , 

6 r 



3.2. LORENTZ STRUCTURE OF THE FUNNEL POTENTIAL 
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results in 



P = 5 1 — 

2 as (r^^) — : a(r ^ 



4 

which implies 

4 

if the Coulomb part dominates, and 

p>2, 

if the linear part dominates. 

Consequently, we arrive at the conclusion that the funnel potential 
V^Yi in Eq. must be a linear combination of a vector and a scalar 
part, 

VMr) = VY{r) + Vs{r) , 
where the Coulomb part Fexch is of vector type, 

Vyir) = V;xch(^) = — , 

6 r 

and the linear part V^conf is of scalar type, 

Vs{r) = Vcon{{r) = ar . 

In summary, from the analysis of the most general conceivable spin 
structures and the experimentally observed quarkonium mass spectra, 
we have been able to determine unambiguously the basic shape of the 
potential acting between quarks: 



The interquark potential Fnr(^) = W(^) + ^s(^) essentially 
consists of a Coulomb part V^xch? which is of vector type, 

6 r 

as well as of a linear part Vconii which is of scalar type, 

Vs{r) = Vconi{r) = ar . 



In this form, the funnel potential represents the genuine prototype 
of all "QCD-inspired" potential models proposed for the description of 
hadrons as bound states of ( "constituent" ) quarks [tl|, 0, ^ • A selection 
of more sophisticated potential models may be found in Appendix |G|. 



Appendix A 



S Matrix, Cross-section, and Decay 
Width 

The normalization of creation and annihilation operators is reflected by 
the (anti-) commutation relations of these operators: 

• For the case of bosons, the nonvanishing commutators are 

[a(p),at(q)]=j(3)(p_q) . 

• For the case of fermions, the nonvanishing anticommutators are 

{6(p, a), 6t(q, r)} = {^i(p, a), ^it(q, r)} = 6^'\^ - q) 5,r ■ 

Normalizing the vacuum state |0) according to 

(0|0) = 1 , 

the normalizations of the one-particle states read 

• for bosons, generically denoted by B, 

(B(p)|B(q)) = j(3)(p-q) 

and 

• for fermions, generically denoted by F, 

(F(p,a)|F(q,r))=(5(3)(p-q)^,, . 

Let us deflne, for any transition i — ^ f from some initial state i to 
some final state f , like, for instance, a scattering or decay process, the 
S-matrix element Sn by 

5fi= (f,out|i,in) = (f|5|i> 
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and the reduced T-matrix element by 

5fi=: Jfi + z(27r)4(5(4)(Pf-Pi)Tfi, 

where Pj and Pf denote the total momenta of initial and final state, 
respectively. The corresponding transition probability Wa is the square 
of the modulus of the transition amplitude Sa — 6{[- 

\Sii - (5fip 



For a finite spatial volume V and a finite time interval T, the obscure 
square of the 6 function may be replaced by 

[{27r)U^^\P{ - POf = {27r)U^^\Pi-Pi}Jd^xexp[iiPf-P0x] 

= {2T^f5^^\Pi-P,)VT, 

which leads to 

Wii = {27r)U^^\Pi - Pi)VT\Tnf . 
The transition rate Rn is the transition probability per unit time: 



Rn^^ = (27r)^ (5(4) (Pf - Pi) V |Tnp 



The cross-section aci is the above transition rate Ra divided by the 
product of the observed flux of the incoming particles, j = n'Ureb times 
the number of target particles, nV; here, 'Urei is the relative velocity of 
the scattered particles and n denotes generically the particle densities. 
If necessary, one has to sum over the final states and to average over 
the initial states, which will be indicated below by a primed sum over 
the possible spin polarizations a: 



1 



an - 



Vrd (nz=l,2 n^) V 

With the particle density 



1 



n = 



corresponding to a normalization volume of size (27r)^, the resulting 
cross-section cr(i f) for the scattering process i f reads 
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APPENDIX A. S MATRIX, CROSS-SECTION, AND DECAY WIDTH 



The product of the energies Ei and E2 of the two scattered particles 
and of their relative velocity Vrei forms a Lorentz invariant: 

El E2 v,ei = \lipiP2Y - m\ ml . 

The decay width Ffi is the transition rate i?fi divided by the number 
of decaying particles nV . If necessary, one again has to sum over the 
final states and to average over the initial states: 




The partial decay width V{i —> f) for the decay of the particle i into a 
particular final state f is therefore given by 

r(i - f) = (27r)7 (^d'p^ 5^'\P, - PO pTnf . 

The total decay width T{i) of the particle i, which is, of course, nothing 
else but the inverse of the average lifetime T[ of this particle, is obtained 
by summing over all possible, kinematically allowed decay channels f: 

r(i)^- = Er(i^f). 



Appendix B 



Feynman Rules for a General 
Gauge Theory 



First of all, a little warning: The correct application of Feynman rules 
requires some experience. In particular, one clearly should be careful 

• when identifying all Feynman diagrams regarded as relevant for 
the specific process under consideration and 

• when computing the combinatorial factors (cf . rule # 3 below) . 

In order to remain on the safe side, it is advisable to evaluate n-point 
Green's functions with the help of "Wick's theorem." Wick's theorem 
allows to convert time-ordered products of field operators, like those 
appearing in the S operator 



5 = Texp ij d^xCi{x) 



into a sum of products of propagators and normal-ordered products of 
field operators. A particular Feynman graph is then nothing else but 
the symbolical representation of a particular operator in the series of 
the Wick decomposition. S-matrix elements may be obtained from the 
Green's functions along the course of the "LSZ reduction technique." 

A general unbroken non-Abelian gauge theory for Dirac fermions '0 
but without scalar bosons is described in gauge by the Lagrangian 

C = -\ Fr + ^iip-m)^-^- {d,V,^f + {d,C) D^C , 

f;. ^ d,v: - + g fabc v: , 

D, = d,-igV^T^ . 

The fermions ip transform according to an arbitrary representation of 
the gauge group. 
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APPENDIX B. FEYNMAN RULES FOR A GENERAL GAUGE THEORY 



The modification of the original theory brought about by the gauge 
fixing must be compensated — in order to maintain the unitarity of the 
S matrix — by adding a further (in general, not Hermitean) term to the 
Lagrangian, which involves anticommuting, scalar "ghost" fields Ca? (a- 

For this theory, the complete list of Feynman rules in momentum 
space for the computation of iTfi, where is the reduced T-matrix 
element defined (in Appendix 0) in terms of the S-matrix element Sn 
by 

Sn = (f, out|i, in) =: 6n + i [2^)^ S^^\Pi - Pj) Tn , 

is: 

1. Propagators: 

Table B.l: Feynman rules for the propagators in a general gauge theory 



Vector-boson propagator 




/i, a ly^l) 




Fermion propagator 


Co = i/j {i p — m) i/j 


i " j 


= i 5'f(/c) , 

Sp(k) — . _ ■ 
f — m + I e K'^ — + z e 


Ghost propagator 




a 5 


= iAF{k)ab , 
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2. Vertices: 

Table B.2: Feynman rules for the vertices in a general gauge theory 



Three-vector-boson vertex 




r, p, c 

1 


= -g fabc [oixy (p - q)p 

+ 9i^p {q - 


Four-vector-boson vertex 


Cl = -\9' fabc fade V^V^'Vj'V: 


\ / 

o 

/\ 
/ \ 

a,d/ \p,c 


^ fi' [/a6e /c(ie {O^J^p 9va 9^(1 g^p) 
"1" face fbde (^Qjiv 9 pa S'/Ucr S'pz^) 
~l~ /orfe /bee (P/U!/ 9a p 9lJ-p 9au)\ 


Vector- boson-rermion vertex 


Cl = gi/jfllj 


a 

1 
1 

^ j 


= i91i.T^i 


Vector-boson-ghost vertex 


Cl = -9 fabc {d.Ca) Cb 


II, c 

1 
1 

p, a ^ ^6 


~ 9 fabcPn 
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3. Symmetry factors: For Feynman graphs involving identical particles 
in internal lines there arise certain combinatorial factors, which have to 
be introduced in order to avoid "double counting." Some examples for 
this are presented in Table pT3 . 

Table B.3: Combinatorial factors for Feynman diagrams involving identical particles 
in internal lines 



Feynman graph 


Statistical factor 


/ \ 


1 


\ / 

\_o_^ 


2! 


/ \ 
\ / 

\ ^ 


1 
2! 


^ \ 

/ \ 




1 
3! 







4. For each closed loop of anticommuting fields, i.e., fermions or ghosts, 
a factor 

(-!)■ 

5. Loop integration: At every vertex, energy-momentum conservation 
has to be taken into account. For every internal and independent four- 
momentum /c, i.e., one which is not constrained by energy-momentum 
conservation at the vertices, an integration 

' (2^ ■ 
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6. External particles: Let the polarization vectors e^(]7, A) describing 
(massless) vector bosons V/j, of four-momentum p and spin polarization 
A be normalized according to 

^e^(p,A) e*(p,A) = -g,,^ . 

X 

Let the Dirac spinors u{p, a) and v{p, a) describing fermions ip of mass 
m, four-momentum and spin polarization a be normalized according 
to 

u{p, a) u{p, r) = Sar , 
v{p,a)v{p,r) = -6ar , 

which is equivalent to 



u\p, a) u{p, r) = v\p, a) v{p, r) = — S^t , 

m 



where 



Ep - + m2 . 

The above normalization implies for the energy projection operators 

1^ + m 



J2u{p,cr)u{p,a) = 

a 

Y.v{p,cr)v{p,a) = 



2 m 
'J) — m 
2m 



With the above conventions, the expansions of the corresponding field 
operators in terms of plane waves read 

• for Hermitean vector bosons 

1 r d^p 



(27r)3/2 



2E. 



J2 A) e^(p, A) exp(-ipx) 



+ a^{p, A) e*(i?, A) ex.p{ipx)] 



'P A 



and 

• for Dirac fermions 



d'^p 



\ 



— [KP: ^) u{p, cr) exp{-ipx) 

rjp cr 

+ d\p^ a) v{p^ cr) exp{ipx)^ 



iIj{x) = 



(27r)3/2 



d^p 



m 



+ b\p, cr) u{p, cr) exp{ipx)^ 
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Table B.4: Feynman rules for the external particles in a general gauge theory 



Particle 


incoming 


outgoing 


Vector boson 








(27r)3/2 pEp"-^^^'^^ 


Fermion 
Antifermion 


1 

(27r)3/2 \ 

1 

(27r)3/2 ^ 


m 

— u{p, a) 
tjp 

m _ 

-Erv{p, a) 
hp 


1 

(27r)3/2 \ 

1 

(27r)3/2 \ 


m _ 

— u{p, a) 
hp 

m 

-Erv{p,a) 
hp 



7. For each change of the relative order of external fermions, a factor 

(-1). 

For illustrative purposes, let's consider electron-positron scattering 
in lowest non-trivial order of the perturbative, that is, loop, expansion 
(cf. Fig. ^): 

e'{pi,ai) + e^{p2,(72) e"(gi,ri) + e+(g2,r2) . 

Applying the above Feynman rules, we may immediately write down 
the corresponding scattering amplitude Tfi (in the so-called Feynman 
gauge, defined by fixing the gauge parameter ^ to the value ^ = 1): 



1 



m 



fi 



(27r)6 ^Ep^ Ep^ E,^ E,^ 
1 



X 



1 



u{qi, Ti) 7^ u{pi, ai) v{p2, (J2) Yv{q2, T2) 
u{qi, Ti) 7^ v{q2, T2) v{p2, CF2) Y u{pi, ai] 
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Figure B.l: Electron-positron scattering in lowest order of the perturbation expansion: 
(a) one-photon exchange, (b) pair annihilation. 



Appendix C 

Quantum Chromodynamics 



Quantum chromodynamics (QCD) is that quantum field theory which 
is generally believed to describe the strong interactions. It is the special 
case of a general gauge theory, characterized by the following features: 

• The gauge group is SU(3)c (where C stands for colour), describing 
the (unbroken) symmetry acting on the colour degrees of freedom. 
The order or dimension of the Lie group SU(A^) is N'^ — 1, which 
equals 8 in the case of SU(3). Hence, SU(3) has eight generators 
T", a = 1, 2, . . . , 8. Of course, the totally antisymmetric structure 
constants fabc^ a,b,c = 1, 2, . . . , 8, are those of SU(3); their values 
may be found in Appendix p.i| . 

• The particle content of QCD comprises the following vector-boson 
and fermion fields: 

— Vector bosons: There are eight, of course, massless, gluons 
transforming according to the adjoint, i.e., eight-dimensional, 
representation of SU(3). 

— Fermions: There are at least six quarks = u, d, s, c, b, t, . . . , 
each of them transforming according to the fundamental, i.e., 
three-dimensional, representation of SU(3). The generators of 
SU(3) in the fundamental representation are given by 

rria _ ^ 

^ fund 2 ' 

where Xa, a = 1, 2, . . . , 8, label the eight Gell-Mann matrices; 
an explicit representation of the latter matrices may be found 
in Appendix |D.4 The total number of quark fiavours will be 
denoted by rip: / = 1, 2, . . . , np. 
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Consequently, the Lagrangian defining QCD reads 

1 flF 

+ gauge-fixing terms 
+ gliost terms , 

witii tiie gluon field strength 

Fflu = ^IJi^u ~ ^v^'li + 9a fabc 

and the gauge-covariant derivative acting on the quark fields 

D, = d,-ig,G;^. 

The parameters of this theory are the strong coupling constant Qs and 
the (current) quark masses rrif. 



Appendix D 
SU(3) 



D.l Structure constants 

The Lie group SU(3) is defined by the following commutation relations 
for its eight generators T", a = 1, . . . , 8: 

[T-,T']=ifabcT'' , 

where the nonvanishing elements among the structure constants fabc are 
listed in Table fn\ . 

Table D.l: Nonvanishing structure constants fabc of SU(3) 



a 


b 


c 


fabc 


1 


2 


3 


1 


1 


4 


7 


1/2 


1 


5 


6 


-1/2 


2 


4 


6 


1/2 


2 


5 


7 


1/2 


3 


4 


5 


1/2 


3 


6 


7 


-1/2 


4 


5 


8 


V3/2 


6 


7 


8 


V3/2 



In addition, the generators satisfy the anticommutation relations 



{T'',T'} = -5ab + dabcT\ 

where the nonvanishing elements among the coefficients dabc are listed 
in Table 0721. 



D.2. GELL-MANN MATRICES 
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Table D.2: Nonvanishing coefficients dabc of SU(3) 





u 
u 


c 




rt 
LL 




C 


f1 

^abc 


1 


1 


8 


l/\/3 


3 


5 


5 


1/2 


1 


4 


6 


1/2 


3 


6 


6 


-1/2 


1 


5 


7 


1/2 


3 


7 


7 


-1/2 


2 


2 


8 


l/v/3 


4 


4 


8 


-l/(2v^) 


2 


4 


7 


-1/2 


5 


5 


8 


-l/(2v^) 


2 


5 


6 


1/2 


6 


6 


8 


-1/(2^3) 


3 


3 


8 


l/\/3 


7 


7 


8 


-l/(2v^) 


3 


4 


4 


1/2 


8 


8 


8 


-l/^/3 



The structure constants fahc are totally antisymmetric whereas the 
coefficients dahc are totally symmetric under permutations of indices. 



D.2 Gell-Mann matrices 



In the fundamental, i.e., three-dimensional, representation of SU(3), the 
eight generators T", a = 1, . . . , 8, are explicitly given by 

fund 2 ' 

where are the eight Gell-Mann matrices 





^0 


1 


0^ 






—i 


0^ 




(I 





0\ 


Ai = 


1 








, A2 = 


i 













-1 







lo 





oj 




lo 





oj 




lo 








/ 1 \ 



1 



A. = 



/ -i \ 




V 



i 



Afi = 



/ \ 
1 



A7 = 



/ 

y i 



^ 
0\ 



-I 




A« = 



1 



/ 1 ^ 

1 
^0 0-2^ 
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D.3 Traces 

The traces of the simplest products of Gell-Mann matrices read 

Tr(Aa) = , 
TTiXaXb) = 2 6ab, ETr(A2) = 16, 

a 

Tr(A« [A^, Ac]) = 4:ifabc , 

TT{Xa{Xb, Xc})=4:dabc , 
Tr(Aa Xb Xc) = 2i fabc + 2 dabc ■ 
Some further useful relations are: 

E {fabc? = 24 , 

^2 40 

{dabc) = , 
a,b,c ^ 

j,k,a ^ ^ ^ 



Appendix E 



We would like to express the second derivatives ViVj$(r) with respect 
to Cartesian coordinates x = {xi,X2,Xs} of an arbitrary function $(r) 
which depends merely on the radial coordinate r = |x| in terms of the 
"spherically symmetric" derivatives coming into question, that is, the 
first and second derivatives of $ with respect to r — which we indicate 
by prime(s) — as well as the Laplacian A = V • V of In other words, 
we would like to rewrite these second derivatives ViVj<l»(r) in terms of 
$'(r), $"(r), and A$(r). To this end, we start from the most general 
ansatz for the expression we are looking for, viz., from 

V,V,*(r) = (aS„ + b^y-<S>'(r) 

where, for every term, the powers of the radial coordinate r have been 
chosen in such a way that the coefficients a, 6, . . . , / are dimensionless. 

It is a simple and straightforward task to determine the coefficients 
a, 6, ...,/: 

• On the one hand, we contract the above ansatz by multiplying it 
by 6ij and by summing over i and j. Using 

5ij 6ij = 6ij 6ji = Tr(l3x3) = 3 , 

we obtain 

A$(r) = (3 a + 6) - $'(r) + (3 c + d) $"(r) + (3 e + /) A$(r) . 
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APPENDIX E. ViVj$(r) 



By comparing both sides of this equation, we may conclude that 
the coefficients a,b, . . . , f have to satisfy the relationships 

3a + 6 = , 
3c + d = , 
3e + / = 1, 



which, in turn, imply 

b 

" = -3' 
d 

^ = -3' 

1 / 
' = 3-3- 

Consequently, taking into account these relations and combining 
corresponding terms, our ansatz simplifies to 



V,Vj$(r) = 



^ J 



b-^'(r) + d^"(r) + fA^(r) 
r 



1 



+ 3%A$(r) 



• On the other hand, for the case i j, we may easily calculate any 
second derivative ViVj$(r) explicitly. The first derivative Vi$(r) 
of $(r) with respect to any of the Cartesian coordinates Xi reads 



Vi$(r) = - $'(r) 



Consequently, for i j, the second derivatives VjVj$(r) of $(r) 
are given by 



ViVj-$(r) = 



$"(r)--$'(r) 
r 



for i^j, 



whereas, by gaining advantage from the fact that 6ij = for i 7^ j, 
the above, already simplified expression for ViVj$(r) reduces to 



V,Vj4)(r) = 



6-$'(r) +£i$"(r) + /A$(r) 
r 



for i j 
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The comparison of these two expressions allows us to fix the three, 
until now indeterminate coefficients b, d, and /: 



6 = -1 , 
d = -hi , 
/ = 0. 



In summary, upon collecting all our findings, the second derivatives 
ViVj$(r) may be expressed as the following linear combinations of the 
"spherically symmetric" derivatives $'(r), $"(r), and A$(r): 



/y* . . I 

" 3 '\ 



$"(r)--$'(r) 
r 



+ i%A$(r) 



Appendix F 

Some Further Formulae for 
Spectroscopy 



For unequal masses mi 777-2, the various spin-dependent relativistic 
corrections to a nonrelativistic potential 

WO = W(r) + Vs(r) 

of vector-plus-scalar Lorentz structure read 

• for the spin-orbit term 

1 1 



^LS = 



4 ml 1712 ^ 



((?77i + 777,2)^ + 2 777,1 1^2) L • S+ 
d 



+ {ml-ml)L-S-]—Vy{r) 

[{ml + m^) L . S+ + {ml - 777?) L • S_] ^^s(r)| , 



with 
and 

• for the spin-spin term 

Hqq = 



S-i- = 81 + 82 
8_ = 81 — 82 ; 



3 mi 7772 



8i-82AW(r) ; 



and 

• for the tensor term 
Ht = 



1 



12?77i 7772 



s 



12 



Id d'^ 

r dr^^^ ^ dr'^^^^ ^ 
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The signatures — parity P, charge conjugation C, and G parity — for a 
quark-antiquark bound state with relative orbital angular momentum 
£, spin S, and isospin / are given by 



P{qq) = 




C{qq) = 


(-1)^+^ , 


G{qq) = 





For instance, for the pion we have £ = 0, 5 = 0, / = ! and therefore 

P(7r) = -1, 
C(7r) = +1, 
G(7r) = -1 . 



Appendix G 

Various Potential Models 



Eichten et al. [] 


1 




4 Qfs 

»/nr = \-ar 

3 r 


Quigg-Rosner 


11 


yNR = ^ln(r/ro) 


Richardson [|T^ 






4 487r2 1 r exp(iq-x) 
3 33 - 2nF(27r)3i^q2 1n(l + q7A2) 


Ono-Schoberl 


13 




{ -- — + d for r < /?i 
V^NR = -6exp(-r/c) + <^ 3 r 

(_ ar for r > Ri 

V 3a Atttq r 
Vs = ar , Vy = I^jr - Vs 


Martin 14 






V^n = A + Br^-' 


Buchmiiller et al. 


|15| 


4 1 /■ ,o , , 47ras(q^) 


Falkensteiner et al. 


4 a 

H^R = -- — erf(A/vrAr) + ar 
3 r 
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Schoberl 17 



4 a 

-erfiAr) + d for r < Ri 

^NR = -5exp(-r/c) + <j 3 r 

ar for r > Ri 



Vs = a r , Vv = Vnr - Vs 



Flamm et al. [IS 



V 



3 (r + ro) 1-107 
4 a, 



3 (r + ro) 1-107 

^NR = W + 



(l-c)+arO-oi(l-rf) 
c + ar°-9irf 
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